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Abstract 

In the context of holographic QCD we analyze Sakai-Sugimoto's chiral model at 
finite baryon density and zero temperature. The baryon number density is introduced 
through compact D4 wrapping S at the tip of D8-D8. Each baryon acts as a chi- 
ral point-like source distributed uniformly over M 3 , and leads a non-vanishing U(V)y 
potential on the brane. For fixed baryon charge density ng we analyze the energy 
density and pressure using the canonical formalism. The baryonic matter with point 
like sources is always in the spontaneously broken phase of chiral symmetry, whatever 
the density. The point-like nature of the sources and large N c cause the matter to be 
repulsive as all baryon interactions are omega mediated. Through the induced DBI 
action on D8-D8, we study the effects of the fixed baryon charge density n# on the 
pion and vector meson masses and couplings. Issues related to vector dominance in 
matter in the context of holographic QCD are also discussed. 



1 Introduction 



Dense hadronic matter is of interest to a number of fundamental problems that range from 
nuclear physics to astrophysics. QCD at finite baryon density is notoriously difficult: (1) 
the introduction of a chemical potential causes most lattice simulations to be numerically 
noisy owing to the sign problem; (2) the baryon-baryon interaction is strong making most 
effective approaches limited to subnuclear matter densities. 

In the limit of a large number of colors N c , QCD is an effective theory of solely mesons 
where baryons appear as chiral skyrmions. Dense matter in large N c is a skyrmion crys- 
tal with spontaneous breaking of chiral symmetry at low density, and restored or stripped 
(Overhauser) chiral symmetry at high density. While some of these aspects can be stud- 
ied qualitatively using large N c motivated chiral models [1], they still lack a first principle 
understanding. 

The AdS / CFT approach has provided a framework for discussing large N c gauge theories 
at strong coupling A = g 2 N c from first principles [2]. A particularly interesting AdS/CFT 
construction is the holographic and chiral approach proposed by Sakai and Sugimoto [3, 4] 
(SS model). In the limit where Nf <C N c , chiral QCD is obtained as a gravity dual to Nf 
D8-D8 embedded into a DA background in 10 dimensions where supersymmetry is broken 
by the Kaluza-Klein (KK) mechanism. The KK scale plays the role of the chiral scale. The 
SS model yields a first principle effective theory of pions, vectors, axials and baryons that is 
in good agreement with experiment [3, 4, 5]. The SS model at finite temperature has been 
studied in [6] and the bayrons in the context of Skyrmion and Instanton have been worked 
out [5, 7]. 

Recently we have suggested that the SS model can be used to analyze dense hadronic 
matter at large N c and strong coupling A [8]. At zero temperature, the quark/baryon chem- 
ical potential \i (or \i% = rriB + N c /i) is introduced as the boundary value of the U(l)v 
brane potential A . The diagonalization of the vector modes in dense matter, enforces vec- 
tor dominance and yields A = fiifj v throughout where ip v is given in [4]. While the mode 
decomposition of ip v on the brane leads to a highly oscillating Aq, we have argued in [8] that 
only those modes in A below the KK scale should be retained. As a result both brane and 
meson properties in holographic and dense baryonic matter were discussed. 

Soon after the posting of this work, several studies appeared addressing the same issue 
of baryonic matter in holographic QCD including also temperature. In [9] it was suggested 
that the Aq field is instead fixed by the equation of motion on the brane by varying the 
pertinent DBI action. In the absence of brane "charges" the authors in [9] concluded that 
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only a constant Aq is a solution, with no baryonic effect at zero temperature. However, the 
conserved baryonic charge has to be mirrored by the "charge" of baryon vertex. The latter 
is obtained by considering brane wrapping of S* 4 within D8-D8. The wrapping number is the 
bulk conserved" charge" which is at the origin of a non-constant A in bulk. 

This point was further developed in [10] and used to discuss the phase structure of dense 
and hot holographic matter, albeit in a brane set up without chiral symmetry. In [11] it was 
argued that the additional "charges" in bulk upset the smoothness of the DBI surface in bulk, 
leading to a spiky structure due to the force balancing condition 1 . As a result, the embedded 
branes should always touch the horizon even for arbitrarily small temperature and/or density, 
thereby altering totally the phase diagram in [10]. This latter point is physically unintuitive. 
Indeed, the "charge" or baryon vertex exists even at zero temperature with no need to 
connect to any horizon. At finite temperature through the insertion of a black hole this 
point is developed in details in [12]. 

In this paper we follow on the analysis in [10] in bulk and at zero temperature (and low 
temperature before deconfinement phase transition). The dual "charge" or baryon vertex is 
inserted at the tip of the minimally embedded D8-D8 surface. This way each of D8, D8 is 
shared equally, leading to a chiral baryon vertex. Since no connecting string is involved, there 
is no spiky structure involved here. Also, there is a one-to-one correspondence between the 
baryon vertex normalization and the Wess-Zumino term in the induced chiral DBI action. 
Thus the boundary chiral skyrmions constructed from the induced DBI action, are dual to 
static and point-like instantons in bulk. These conditions will be relaxed in a sequel. In 
many ways, this approach complements the original discussion in [8]. 

In section 2, we briefly review the SS model and set up the notations. In section 3, we 
introduce the U{l)y field A in bulk and show how the baryon charge density n B affects its 
minimal profile. In section 3 and 4, we construct the bulk hamiltonian and derive the energy 
density as a function of the identified baryon density. The energy density is found to grow 
about quadratically with the baryon density. In section 5, we summarize the construction of 
the chiral effective action for pions, vectors and axials at zero density. In section 6, we show 
how this chiral effective action is modified by the finite "charges" in bulk. A number of meson 
properties are discussed as a function of the identified baryon number. Our conclusions are 
in section 7. Throughout, the canonical formalism will be used. 

1 Notice that while the baryonic charges are space separated, they occupy the same position in the dual 
transverse space! 
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2 SS Model 



In this section we summarize the D4/D8-D8 set up for notation and completeness. For a 
thorough presentation we refer to [3] and references therein. The metric, dilaton 0, and the 
3-form RR field C 3 in N c D4-branes background are given by 

ds 2 = (^y 2 (v^dx» + f(u)dr>) + (^) 3/2 ^ + uHn\ 

e^(£) 3/ \ F 4 = dC 3 = 2 -^e 4 , m = l-?k,, (1) 

where x M = x 0,1 ' 2,3 , r(= x 4 ) is the compact variable on S 1 . U(> Ukk) and f2 4 are the radial 
coordinate and four angle variables in the x 5 ' 6 ' 7 ' 8 ' 9 direction. R? = rcg s N c l 3 , where g s and l s 
are the string coupling and length respectively. V4 = 87r 2 /3 is the volume of unit S 4 and e 4 
is the corresponding volume form. 

To avoid a conical singularity at U — Ukk the period of 5r of the compactified r direction 
is set to 

^KK 

in terms of which we define the Kaluza-Klein mass as 

_2vr_3< 2 

Mkk = 57 " 2^ • (3) 

The parameters R, Ukk, and g s may be expressed in terms of M KK , A(= gyMN c ), and Z s as 

B ' = 5i^' ^K^AM KKi ;, k-^jj^ (4) 

Now, consider Nf probe D8-branes in the N c D4-branes background. With U(Nf) gauge field 
Am on the D8-branes, the effective action consists of the DBI action and the Chern-Simons 
action: 

o qDBI , Q cs 

>->D8 — >JD8 + D D8 > 

SE 8 BI = -T 8 J d 9 x e-* tr y/-det(g M N + ^a'F MN ) , (5) 
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where T 8 = l/((27r) 8 /;?), the tension of the D8-brane, F M n = <9m A n — 9nA m — % \A M -,A N \ 
(M, N — 0, 1, • • • ,8), and g M N is the induced metric on D8-branes: 



ds 2 = ( - J ri^dafdx" + 



^) 3/ w)(T'(u)r+(^ 32 



r) ww>r + [jj) W) 



R 



3/2 



dU 2 + ( jj ) U 2 d£l\ , (7) 



where r' — The effective action of D8 has the same form and the total action of Nf 
D8-D8-branes has a symmetry 

U(N f ) L x U(N f ) R = SU(N f ) L x SU{N f ) R x U(l) v x U(1) A , (8) 

which is interpreted as a flavor chiral symmetry of massless quarks. 

In the SS model baryons are skyrmions in M 3 . However, in the gravity dual they are 
either an effective fermion degree of freedom in 5-dimensions (bottom-up) or an instanton 
wrapping DA (top-down) and sourcing the baryon current through the Chern-Simons term 
[5, 8, 13, 14]. In both cases, the baryon can be treated as a delta function source of the 
brane gauge field A (which is u in M 3 ), which we will use explicitly. 



3 Background field A 



o 



Let Aq{U) be a U(l)v valued background gauge field in bulk. Its boundary value is related 
to the baryon chemical potential [8, 9, 10, 12]. In the absence of the source, the effective 
action of the D8-branes (5) becomes 



5e>8 — 



NfT * V < I d'xdUU* 



/(r ' )2+ (f) (^-Wf) 



(9) 



where A' = ^f- and the Chern-Simons action vanishes. The equations of motion for t(U) 
and Aq(U) are [9] 



d 

dU 



U 4 fT> 



0, 



(10) 
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d 

dU 



_SJ 7 iff + (§Y (f- 1 ~ (2na'A' ) 2 ) 



(11) 



In this paper we consider only the case r' = 0, Sakai-Sugimoto's original embedding [3, 4], 
where the D8-branes configuration in the r coordinate is not affected by the existence of 
background Aq. This corresponds to r = j, the maximal asymptotic separation between 
D8 and D8 branes. 

To compare with [3, 4] we change the variable U to z through 



U= (U^ K + U KK z 2 ) l/3 



(12) 



The action (9) is then 



S 



D8 



-N f T I d 4 x j™dzU 2 ^ 1 - (2irarl^(d z A y , 



(13) 



where we used r' = and T = 2 ^^^ 3 . It is useful to define the dimensionless quantities 



KK 



-, K(U) = 1 + Z 2 



U 



E/i 



KK 



(14) 



in terms of which the action is written as 2 



Sds — —a 



j d A x J dZK 2 ' 3 ^l-bK^idzAo) 2 , 



where 



3% 5 



3 b 7T 



4A 2 M 2 K 



(15) 



(16) 



Now we introduce the baryon source coupled to Aq through the Chern-Simons term 
[8, 5, 13] as mentioned before. We assume that baryons are uniformly distributed over R 3 
space whose volume is V. For large A, the instanton size is l/y/X [5, 14]. It can be treated 
as a static delta function source at large iV c . For a uniform baryon distribution, the source 
is 

Source = N c n B J d A x J dZ5(Z)A (Z). (17) 



2 The integral is extended to (— oo, oo) to take into account D8 branes as well as D8 branes. 
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The equation of motion of Aq is 

d dL 



dZd(d z A^ '' 

which yields 

dL 1 



n q S(Z) , (18) 



d(d z A ) 2 



n q sgn(Z) , (19) 



where n q = N c tib is the quark density and the step function sig(Z) is determined by the 
symmetry between D8 [Z > 0) and D8(Z < 0). By integrating once more we get the classical 
solution A Q 

A (Z; n q ) = A)(0) + / dZ - ql . (20) 

Jo J(ab) 2 K 2 + bK 1 / 3 n 2 q /4: 



We introduce the "baryon charge chemical potential of a quark", /i, by [10, 12] 

fj,(n q ) = lim A {Z;n q ) . (21) 

\Z\— >oo 

This relation also defines /i as a function of n g and vice versa. Furthermore we define the 
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baryon chemical potential as 



fi B = m B + N c fx . (22) 

In Fig. (la) we plot the profile of A (Z) in the Z coordinate and in Fig.(lb) we show /j, for 
various baryon densities. Since we work in the canonical formalism u is more like a Lagrange 
constraint. 

Throughout, the numerics will be carried using the following values [3, 4]: Nf — 2, 
N c — 3, fir — 92.6MeV, and m p = 776MeV. The smallest eigenvalue was calculated to be 
Ai = 0.669. Using these five values we can estimate Mrk, A, k, a, and b: 

m 547r 4 XN 

M KK = -^=^ 950MeV, A = g 2 YM N c = fljfj^- - 16-71, « = ^ 0.0075, (23) 

and 

a = 3.76 • 10 9 MeV\ b = 7.16 x lO^MeV^ 2 . (24) 

The definition of k and A are different from [3, 4] by a factor of 2, but it is consistent with 
[5]. In all figures is normalized to with no the nuclear matter density, 

n = 0.17/m" 3 ~ 1.3 x !0 6 MeV 3 . (25) 



4 Thermodynamics 

Consider the action (15) with the source term (18), 

/r+oo 
d A x I dZL 
J —oo 

with £j = -aK 2 / 3 l-bK l /\d z A ) 2 + n q 5{Z)A {Z) . (26) 

The Aq is an auxiliary field with no time- dependence. It can be eliminated by the equation 
of motion (19) and (20). The energy is 

+oo 



/r+oo 
dx 3 dZ{-L) 
J — oo 



/+oo / „2 

dZK 2 ' 3 \jl + - n g u , (27) 
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where V is short for f dx 3 and we may set /i — 0. The chemical potential ji is constrained 
by the Gibbs relation \x = 9F q^ where F(n q ) is the Helmholtz free energy which is U (n q ) 
at zero temperature. Thus 



-L 



dZ 



n q /A 



oo ^J(ab) 2 K 2 + &ATi/3 n 2/ 4 



(28) 



which is in agreement with the solution (20) for Aq(0) = 0. We note that this construction 
is consistent with [8, 9, 11, 10] where the grand potential is identified with the DBI action 
at finite /i. 

In terms of the baryon number density n B (n q /N c ) the regularized Helmholtz free energy 

is 



Fre g (n B ) 
aV 



- / 

J — c 



dZK 2 ' 3 



1 + ( N * n ?) 2 K -W - 1 



Aa 2 b 



(29) 



after subtracting the vacuum value. The regularized internal energy U, pressure p and grand 
potential Q as a function of baryon number density n B or the baryon chemical potential \x B 
are 



U reg (n B ) 



aV 



p(n B ) 



reg 



^reg^fi) 
Ol 7 



fJ-B 



/oo 
-oo 

30 

dZK 2 / 3 

oo 

30 



i + {N ^ b)2 k^ - 1 



4a 2 6 



1 - 



^ 1 + ^^-5/3_ ' 

1 

= -1 

/ j _|_ (Af e rt B (^ B )) 2 ^_5/3 

-oo y/(ab) 2 K 2 + bK 1 /3(N c n B /2) 2 



(30) 



where JTs = fi B — m B = N c /i. 

In Fig. (2) we present the numerical plots of these thermodynamic functions with the 
numerical inputs in section 3. For small baryon densities the energy density is quadratic in 
n B /no (or fi/m p ). At large baryon densities it is of order (n B /no) 1A - The small density limit 
can be qualitatively understood by noting that in bulk the A configuration for fixed charge is 
obtained by minimizing the induced DBI action of D8-D8. Thus only flavor-meson mediated 
interactions between the point-like baryons are included. At large N c the DA mediated 
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Internal energy density 



Pressure 




Figure 2: Numerical behaviour of the thermodynamic functions: See Eq.(30) 



correlated gravitons (glueballs on the boundary) are heavy and decouple. Since our point 
baryonic vertices in bulk map on infinite size skyrmions at the boundary this implies that only 
u} exchanges survive at large N c . Rho and pion exchange relies on skyrmion gradients which 
are zero. At low baryon densities, the dominant Skyrmion-omega-Skyrmion interaction is 
two-body and repulsive. Thus the energy density is positive and quadratic in the baryon 
density. The baryonic matter is prevented from flying apart by the container V. At large 
baryon densities, the energy density softens as the quark chemical potential is seen to saturate 
to (n B /n ) 0A numerically. We recall that the baryons are fixed sources so no Fermi motion is 



Thermodynamic function 


n B /n ~ 


n B /n ~ 10 


n B /n — > oo 


Internal energy 
Pressure 
Chemical potential 
Grand potential 


(n B /n ) 2 
(n B /n ) 2 
(n B /n o y 
-(jr B /m p ) 2 


(tib/tio) 1 - 85 
(n B /n ) 1A5 
(n B /n )°- e7 
-(jr B /m p ) 2M 


(n B /n ) 1A 
(n B /n ) 1A 
(n B /n ) 0A 



Table 1: Numerical behaviour of the thermodynamic functions: See Eq.(30) 



9 



involved to this order. The pressure behaves as (n^/no) 2 at low baryon densities, and again 
softens to (n B /n ) 7 ^ at large baryon densities from the plot. We summarize the behaviour 
of the thermodynamic functions obtained numerically in Table. (1). In this paper we do not 
consider the back reaction of gravity for baryons or D8 brane, therefore the behaviour at 
higher densities, say ns/n ^> 10, is not justified. 

5 Effective meson action: ub = 

In [3, 4] the meson spectrum and coupling was studied at zero baryon density by analyzing 
the DBI action of D8-D8 branes with the fluctuating gauge field Am- We want to extend the 
analysis to finite baryon density or n B 7^ 0. For this purpose we streamline in this section 
the construction in [3, 4] for notational purposes and completeness. In the next two sections 
we add the background U{l)y field A to the fluctuating gauge field A M . It will enable us 
to study meson properties at finite baryon density. 

5.1 Mode decomposition of A M 

The gauge field Am has nine components, A^ = ^2,3,4, A z {= A 5 ), and A a (a = 5,6,7,8, 
the coordinates on the S* 4 ). We assume that A a = 0, and A^ and A z are independent of the 
coordinate on S* 4 . We further assume that A M can be expanded in terms of complete sets, 
ip n (z) and (f) n (z) as 

00 

A^,z) = (31) 

?1=1 

OO 

A x (x»,z) = ^ ( °V)0o(*) + £> (n V)0n(z) , (32) 

71=1 

where is identified with vector and axial vector mesons and (p^ with pions. ip^ can be 
absorbed into through the gauge transformation (section 5.2). ip n satisfies the eigenvalue 
equation, 

-K 1 / 3 d z (Kd z ^ n ) = X n ^ n , (33) 

with the boundary condition (9^^ n (0) = (vector meson) or if) n (0) = (axial vector meson) 
at Z = 0. They are normalized by 

k I dZ K-V 3 M m = 5 nm , (34) 



10 



where re = T(27ra') 2 R 3 = ^gfe , and (33) and (34) implies 

k J dZ K(d z 4) n ){dzi>m) = KA im • (35) 
The 4> n (Z) are chosen such that 

UZ) = * d z MZ) (n>l), (36) 

VA n IVl KK UKK 

MZ) = jntou/r KK Tc ' (37) 



with the normalization condition: 

(0m, <Pn) = ^Ik^KK J dZ K 'Pm'Pn = &mn , (38) 

which is compatible with (35). 

5.2 Effective meson action 

With the gauge field A^x^.z) and A z (x^,z) the DBI action of the D8-D8-branes becomes 
5-dimensional 3 : 

= -f J d 4 xdz U 2 

I R 3 9 U 

tr J 1 + (2na') 2 -—F lxv F^ + (27ra'y--—F flz F^ + [F3] + [F 4 ] + [F 5 ] ,(39) 

y 4 (7kk 

where T = ^^F, U is a function of z by (12), and the indices are contracted by the metric 
(— , +, +, +, +). [F 3 ], [F A ], and [F 5 ] are short for the terms of F 3 , F A , and F 5 respectively 
Notice that the range of z is extended from [0, oo] to [— oo, oo] to account for both D8 and 
D8. 

Inserting (31) and (32) into (39) and using the orthonomality of ip n and 4> n ((34)~(38)), 
we have [3, 4] 

SSSs ~ [ *xtr \(d^) 2 + f; (\{d,B^ - d v B^) 2 + X n M 2 KK (B^ - X^d^f 

J n=l ^ 



+ (interaction terms) . (40) 



3 The gauge group generators t a are normalized as trt a t = 5 a b/2 
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Here and are interpreted as a masseless pion field and an infinite tower of vector 
(or axial) vector meson fields with masses m 2 n (= A n M£ K ). The lightest vector meson p is 
identified with p^ are absorbed into In the expansion (31) and (32), we have 

implicitly assumed that the gauge fields are zero asymptotically, i.e. Am(x^,z) — > as 
z — > ±oo. The residual gauge transformation that does not break this condition is obtained 
by a gauge function g(x tt ,z) that asymptotes a constant g{x^,z) — > g± at z ± oo. (g + ,g_) 
are interpreted as elements of the chiral symmetry group U(Nf) L x U(Nf) R in QCD with 
iV/ massless flavors. 

5.3 A z = gauge and pion effective action 

In the previous subsection we worked in the gauge A M (x^, z) — > as z — > ±oo. However the 
A z = gauge can be achieved by applying the gauge transformation A M — > gA M g~ Xj rgdu9" x 
with the gauge function 

( ? - 1 (^,^)=Pexp|-^^ / ^(^,z / )} . (41) 
Then the asymptotic values of A^(z — > oo) do not vanish and change to 

A^, z) - id^)%t±(??) as ±oo , (42) 
where £±(x M ) = lim 2 ^. ±00 g(x^ ) z). The gauge fields can be expanded as 

oo 

a^,z) = e + (^)^c 1 (^)(^)^ + (^)+e4^)^C 1 (^)^4^) + E 5 i n) ( xAt )^W > 

n=l 

A g (x»,z) = 0, (43) 

where tp± is the non-normalizable zero mode of (33) with the appropriate boundary condition 
to yield (42): 

i>± = o =•= & , 

■00 = - arctan(Z) (44) 

7T 

There is a residual gauge symmetry which maintains A z = 0. It is given by the z- 
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independent gauge transformation h(x >l ) 1 



A M (x», z) -> h{x^)A M (x^ z)hr\x») + h{x»)d M hr\x») , (45) 

which acts on the component fields as 

e± - ht ± g±\ (46) 
_> hBtofT 1 , (47) 

where we considered chiral symmetry (? ± together. Then £±(x M ) are interpreted as the U(Nf) 
valued fields ^l,r(^) which carry the pion degrees of freedom in the hidden local symmetry 
approach . Indeed the transformation property (46) is the same as that for ^r^) if we 
interpret h(x^) G U(Nf) as the hidden local symmetry. They are related to the U(Nf) 
valued pion field U(x^) in the chiral Lagrangian by 

^(x^-ix^ 1 ) = U(a?) = e 2m ^ f " . (48) 

The pion field Il(x^) is identical to (p^(x^) in (32) in leading order. A convenient gauge 
choice is 

Z_( x ») = 1, e; 1 ^) = U(x^) = e 2m{x " )/f * (49) 
which expresses the gauge fields as, 

A^,z) = U-\x»)d,U(x^ + (z) +Y,Bp\x»)M*) (50) 

n>l 

In this gauge, after omitting the vector meson fields the effective action reduces to the 
Skyrme model 



qDBJ_ 

°D8-D8 



,(»)_ 



I d'x (^tr (U-%U) 2 + JLtr [U-^U-^U] 2 ^ , (51) 



where e s 2 = k Jdz K 1 / 3 (1 — -?/> 2 ) 2 and the pion decay constant f n is fixed by the comparison 
with the Skyrme model: 

ft = ^ M kk = ^M 2 KK \N C , (52) 
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Another gauge we will consider below is 



^(x") = C-(^) = e in( * M)// ^ . (53) 
in terms of which the gauge fields are written as 

oo 

A^x^z) = a^oW+^+^fMiW, (54) 

71=1 

«„M = {r 1 ,^} = ^n+[[^n 3 ]] + o(n 4 ), 
= ^[r 1 ,^] = ^[n,^n] + o(n 4 ), 

where [[c^n 3 ]] = ((<yi)n 2 + rr%n - 2n(^n)n). 
6 Effective meson action: ub ^ 

We now extend the previous analysis to finite baryon density for % = 0. This is achieved 
by adding the background U(l) v field A to the fluctuating gauge field A M . Since, the 
vacuum modes {ip n , <p n } are not mass eigenmodes in matter, we may choose more pertinent 
eigenmodes in matter. Two basis set are possible: (1) medium mass eigenmodes ip n ~ 
e~ imt fn(z); (2) screening eigenmodes ip ~ e lk ' x f n (z). With this in mind, we have the following 
gauge fields decomposition 

oo 

A (x»,z) = A,{z) + Y,Bt\x^n{z) , (55) 

rt=l 

oo 

Mx»,z) = Ys B i n) ( x ")Mz) , (56) 

71=1 

OO 

A g (x»,z) = 5> ( "V)0nW. (57) 



n=0 



A (z) is the background gauge field. The time component modes (cu n (z)) and the space 
component (ip n (z)) are not necessarily the same as Lorentz symmetry does not hold in the 
matter rest frame. Note that F^ z is modified by ^Iq while F nu is not. 
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2d Q Bl n) d B ( - m ^ifj n ifj m + 2d i B { ™ ) d j B 
-2d Bi n) d l B^ n ij m 



Ho, 



+{diBf ] - djB\ n) )(&BWi _ gj fiH*)^^ 



>(n) 



= a 2 



F„ Z F^ 



-{A f 



+2A Q 



A) 



+ 



,9oV? (n) <9V m) 0n0 m + B^B^UnOJrn - 2do<fW B™ <j>„ 



+( 9 iV 9( n )9V {m) 0n0 m + B^B^M m - 2 V n) £ (m)i <M>J = (3, 



[F 4 



.iu.r j (M 2 ci 



= 72 



Table 2: The relevant terms in evaluating the DBI action up to quadratic order in the fields 
(B^,ip). The upper dot stands for the derivative with respect to z. The terms should be 
understood in the integral and trace operation. 



In order to compute the DBI action (39), 
-f J d A xdz U 2 



qDBI_ 

D8-D8 



i? 3 9 U 

tr x \ 1 + [2na'f^F, v F^ + (2ira<f-—F, z F^ + [F 3 ] + [F 4 ] + [F 5 ] , 

we need to know F^F**", F^F^, [F 3 ], [F 4 ], and [F 5 ], which are involved in general. To 
quadratic order (ignoring 0((B^, y?) 3 )), the contributions are greatly simplified because of: 
1) cyclic property of the trace, 2) antisymmetry of F^ u , 3) parity of mode functions. Then 
there is no contribution from [F 3 ] and [-F 5 ]. [F A ] has important terms that will modify 
F^: 



IF 4 } = (2^— (f ) 3 F Qz F°*F l3 FV + QUB^Y 



(58) 



Table (2) lists all the relevant terms, where we have introduced defined as 
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fij = diVj - djVi , 



(59) 



with i,j = 1,2,3. Table (2) should be understood in the integral and trace operation. We 
omitted some terms vanishing in those operations and rearranged some terms by using the 
cyclicity of the trace. 

In terms of the definitions on the RHS of the Table (2) , the action reads 



qDBI_ 
°D8-D8 



— T J d 4 xdz U 2 tr y/Po + P! , 



with 



Po = i_(W) 2 H^ft = i-wd(0 z 4>) 2 , 



4CA 



KK 



Pi = (2W) 2 ^(« 2 ) + (2W) 2 ^(A+ft) 



+ (2W) 



,. 4 9 R 3 



8U KK U2 



(72) , 



(60) 



(61) 



(62) 



where P does not contain meson fields but involves the baryon density. Expanding the 
action for small fields we have 



^DBJ_ 
'D8-D8 



— T / (fxdzU 2 tr 



Po + 



1 Pi 1 P 2 



2v^ 8^ 



+ o((p^) 3 ) 



= S 1 + S 2 + 0((B^vf) , 



(63) 



with 



S 1 = -T J d A xdzU 2 tr A' 1 , 
So = — T / r/ ; .rr/: C/ 2 tr 



-APx - ^A 3 P 2 
2 1 8 1 



(64) 
(65) 



where the modification factor A(n#) is 



A(n B ) 



^ yJl-bK*(d z Ao¥ 



1 + 71^- 5/3 
4a 2 6 
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-So is the grand potential discussed in section 4, and 5*2 will be reduced to 4 



So = -tr 



+ 
+ 
+ 



J dZK~ 1/3 A * n tt„ 



d Bl m) d°B^ 1 + 
2d B<; n) d i B im)0 



J dZK~ 1/3 A tt n tt r 



diB^d'B^ 



-{diBf ] - djB\ n) ){j^B^ m)i - d j B (m)i ) 
2 J 



Ml K I dZKA 3 d z n n d z a r 



M| K / dZKA dz^ndz^r 



Ml K I dZKA^A 



I 



Ml K / dZKAH n d z Sl r 



d tp {n) d°ip {m) + 
2d <p {n) B (m)0 



Ml K j dZKA<S> n <S> m 



M£ K / dZKA^dz^, 



2cV n) 5 (m)i j , 
(66) 



where we defined the scaled eigenfounctions as 



(67) 



At zero density A = 1, so <J» n = Q n and the action reduces to the (40) by the same mode 
function in (33) ~ (38). However at finite density the eigen modes fl n , \l/ n , and <3> n cannot 
be determined uniquely. In other words there is no mode decomposition which makes the 
action completly diagonal. So we consider the space-like and time-like separatly: (1) A M = 
A M {x\z) and (2) A M = A M (x°,z). 



4 Note that the pattern: A, A 1 , and A 3 . This pattern appears also when we consider higher order terms 
including couplings. The origin is explained in Appendix A. 
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6.1 Space-like fields Am = Am(x\ z) 

First we consider time-independent gauge fields. Up to quadratic order the action is 



S 2 = —1 v i 'i'r 

+ 
+ 



dZK- 1/3 A n„ a 



M| K J dZKA 3 d z Q n d z Vt 



-{diB^ - djB^i&BW' - d j B (m)i ) 
2 



M^ K / dZKA dz^ndzV 



(n) r>(m)i 



BT'B 



+ 



Mi K / dZKA^ s m 



M| K I dZK A&z^i 



■ n * m 



2dw {n) B i{m) \ , 



where we have defined the scaled eigenfunctions as 



(68) 



To diagonalize the action we choose as the eigenfunction satisfying 

-K^A^dz (ka 3 d z n n ) = \%n n , 

-K^A d z (KA d z * s n ) = XX , 



(69) 



(70) 
(71) 



with the normalization conditions, 



J dZK^ 3 A n n n m = 5 nm , 

J dZK-^A-^ s m = 5 nm , 



which imply 



/ 
/ 



dZKA 3 d z tt n d z tt m 
dZKA dzKdz^l = \ S J r , 



/ ^ n u nm 1 



(72) 
(73) 



(74) 
(75) 
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If we choose <&:? as 



n Mkk \/Aif 



d z * 8 n (n > 1) , 



1 



Mi 



KK 



y/fdZ(K-iA 



(76) 



then di(p( n ^ (n > 1) can be absorbed into -B^™' 1 through the gauge transformation 



1 



Mkkv^I 

These choices of mode functions reduces the action to 



(77) 



S 2 = -tr y d 4 :r jV 0) <9V {0) 

+ ( 9, J Bj n) ^ J BW° + ^4 n) / ( " )ij + K 2 £o°£ (n)0 + M^ 2 B{ n) BW*} , (78) 

where we have defined longitudinal screening masses and transverse screening masses 
Mf as 



6.2 Time-like fields A M = A M (x°, z) 
For spacially homogeneous gauge fields the action reads 



KK 



(79) 



+ 



+ 



/ 



M£ K / dZKA d z W n d z ^ 



+ 



Ml K j dZKA 3 d z n n d z n m 



I 



Ml K I dZKA 3 ^^ 



Ml K J dZKA 3 ^d z n r , 



2d v {n) B°( r ' 



(80) 



where we have defined the scaled eigenfunctions 



n 



(81) 
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We choose as the eigenf unction satisfying 



-K^A^d z (KA d z * T n ) = XX , 
-K^A-'dz (KA 3 d z n n ) = A^ n , 



with the normalization conditions, 



which imply 



J dZK-^A = 5 nm , 

J dZK'^A n n n m = 5 nm , 

dZKA d z *Zd z *l = \ T n d nm , 

J dZKA 3 d z n n o z n m = \% m , 



If we choose <£>:? as 



then d (f^ (n > 1) can be absorbed into Bq through the gauge transformation 



(n) 



yJfdZ(K-iA-*) A ' A:i 



1 



Mkkv/AJ 



The action is reduced to 



S 2 = -tr J d 4 x{d ^ n) d° V {n) + d Bl n) d°B^ + mlB^BW + M^A^^W } , 



where we have defined the mass 



KK 



6.3 Pion effective action 

In this subsection we work in the A z = gauge following the procedure in section 5.3. 
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6.3.1 Time-like field (A M = A M (x°, z)) 

First consider the case Am = Am(x°,z). By the gauge transformation Am 
gdM9~ l with the gauge function 



g-\x\ z)=P exp j- jT dz' A z (x°, z') J , 



gA M g 1 + 



(92) 



the gauge fields are rewritten as 



A (x°,z) = yi W+e + (^ )5oC 1 (^V + W+^(^ )5oC 1 (^ KW , (93) 
Mx°,z) = A z {x°,z)=0, (94) 

where we have omitted the vector mesons The uj± are obtained as zero mode solutions 

of (83) satisfying the boundary condition for A (x°,z): 



2 JdZiK^A- 3 ) J KA 3 ' 



dZ- 



(95) 



By using the residual gauge symmetry h(x^) (49) and (50) we may express the gauge field 
as 



A (x°, z) = A + U- 1 (x G )d U(x' ) )uj + (z) 



(96) 



The field strength is 



F z ^A + U-'d U^(z) , F IU , = 0, 



(97) 



where 



4%{z) = d z u + (z) 



U KK f dZ(K- l A-*) KA 3 ' 



The action becomes 



(98) 



So = tr / d 4 x 



KK 



JdZK^A- 



(u-'douy , 



(99) 



21 



and we identify the time-like pion decay constant fj as 



f 2 

J 7T 



JdZK^A' 3 ' 



(100) 



by comparison with the Skyrme model. 



6.3.2 Space-like field (A M = A M (x\z)) 

Similarly we consider the case A M = A M (x\z). Using the same gauge transformation we 
can work with the gauge fields, 

A( x \z) = e + (^)^; 1 (^)^+W + e4^)^c 1 (^)^W , 

A (x\z) = A (z) , A z (x\z) = , 

where are obtained as a zero mode solution of (73) satisfying the pertinent boundary 
condition of A^x 1 , z): 



JdZ(K-iA-r) I dZ lk 
Then the gauge field and the field strength in the gauge (49) are 

A (x°,z) = Ao + U-'ix^Uix 1 )^) 
F z , = Ao + U^dWKz) , 

where we do not consider F^ v since we are interested in the kinetic part and 



^ (z) = d z ^ + (z) 



U KK $ dZ(K- l /\- r ) KA ' 



The action is 



So = tr / d A x 



kM, 



KK 



jdZK~ l A- 1 



(tr 1 w , 



and f£ is identified by 



J TV 



AkM, 



KK 



JdZK^A- 1 ' 



(101) 



(102) 



(103) 



(104) 



(105) 
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6.4 Vector Mesons Interactions 

In this section we study the interactions of the fields B^\ and corresponding to 
the lowest medium modes Q±, and $1. For simplicity, we use the following notation, 

v = , Vi = , n = v? (0) . (106) 

The details of the computation are relegated to Appendix A 5 . 

6.4.1 Time-like Fields A M = A M (x°,z) 

S 2 = tr J d A x{ - d Ud°U + d v l d°v i + m\viV l + Ml K \°v v } 

-2g T vIl , vo[IL, d°U] + g T v , 2d v i [v°, v 1 } + • • • } , (107) 

where the couplings can be read from (4,5) in Appendix A by substituting the vacuum mode 
functions by the medium mode functions 

<& = 4= I dZK-^n^JfA. (ios) 
V K J 

6.4.2 Space-like Fields A M = A M {x\z) 

S 2 = ti J d A x{ - diUdm + d l v d l v° + ^fijf ij + Mf v v° + M^ViV 1 

-2g s vU2 Vi [U, dm] + g s vi 2d iVo [v\ v°] + g% f l3 \v\ j\ + ■ ■ ■ } , (109) 

where the couplings can be read from (4,5) in Appendix A, again by substituting the vacuum 
mode functions by the medium mode functions 

5 Both in Appendix A and this section, the vector meson field is considered as anti Hermitian. Although 
in Appendix A, we are working with the vacuum modes instead of the medium modes, the conversion can 
be done by inspection using the formula tabulated in Table (4,5) 
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S _ j J KA 

g s v s = -L [ dZK-^^ffA' 1 . (110) 
v K J 

6.4.3 Zero Density Limit 

To check the current mode decomposition used in this section, we take the zero baryon 
density limit. In this case, Lorentz symmetry is enforced and the action reads 

S 2 = tr J d 4 x{ - d^U + \f^r + mWv"} 

-2g vU 2 v^U, d»Ti\ + g v s f^v", v"] + • • • } , (111) 

which is the same as Eqn.(5.40) in [3] except the vUU coupling. The difference comes from 
the gauge choice. In [3] A z = gauge is used and we chose A M (z — > oo) — > 0. Since the 
difference is merely a gauge choice, physics will not be changed. However we will repeat 
the analysis of couplings at zero density with (111), since the action in our gauge is more 
convenient for reading off physical quantities. Also it is readily extendable to finite baryon 
density. 

First we examine the KSRF relation by defining czksrf as 
_ 4 ^/' J 2 - 03 Experiment 

OKSRF = 2 ~ \ ' 

J^i [1.3 Sakai Sugimoto model 

which is the same value reported in [3, 4], as expected. The universality of the vector meson 
coupling can be checked by au defined as 

g vIl 2 I 1 The universality of the vector meson coupling 

au = ~ <^ , (113) 

gv 3 [ 0.93 Sakai Sugimoto model 

which is also the same value as in [3]. Notice that both relations include g vU _2 and can be 
read from (111). In the A z = gauge we should convert o^n 2 to a vT] 2 by [3] 

a«n2 = — • (114) 
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When we consider the field redefinition in (111) 

^-^ + ^[n,^n] , (us) 

the algebraic relation (114) appears immediate. However when we look at the integral 
expression of g V Yp and a„n 2 the equivalence is obscured. 

_ 2g vI p 

a vU 2 — 2 



mi 



s 2 



y J dZK-^^l-A^l) J dZK(dz^i) 2 = J dZK~H x . (116) 

Next and following [3], we compare (111) with the action from the hidden local symmetry 
approach 

S H = tr J d'x[ - ^ncrn + l -f, v r + ag 2 f w v^} 

-agv,[Y[,d^}+gf^y} + ---] . (117) 

The hidden local symmetry parameter (LHS) can be written in terms of the D-brane effective 
action parameter (RHS): 

9 = 9v* , (118) 

a = ^ = ^, (119) 

9 g v z 

ag l 2g v3 g vU 2 

where we used the first two relations to get the last. We may define the parameter an which 
quantify the difference between hidden local symmetry approach and our model: 6 

_2g v sg vIl 2fl fl Hidden local symmetry 

a H = 5 ~ < , (121) 

mf [0.72 Sakai Sugimoto model 

which is the same value reported in [3], as expected, an may be interpreted as follows. Since 
f n is an input parameter the Hidden local symmetry has two adjustable parameters, so a 
is not uniquely determined. It can be fixed by (119) or (120). When these two procedures 



6 For a=2 the hidden local symmetry approach implies KSRF relation and the universality of the vector 
meson coupling. Here, we do not require this value since we want to compare our model with the hidden 
local symmetry itself. 
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yield the same value, a\j — 1. 



6.5 Numerical results 

All the numerical work reported here has been carried out for the lowest modes v 
and n = (pW with the parameters discussed in section 3. 

6.5.1 Mass and Screening Mass 

From the previous section the meson masses (91) (time-like) and the screening masses (79) 
(space-like) are defined as 

m n = \AT m kk , 

= v^Mkk , K = y^M KK , (122) 

where A^, A^, and Af are determined as the eigenvalues of the following equations ((82), (83), (73)), 
respectively: 

-K^A-'dz (KA d z *l) = XX , (123) 
-K^A-'dz (KA 3 d z Q n ) = A^ n , (124) 
—if 1 / 3 A d z (KA d z *%) = XX ■ (125) 

Their dependense on the baryon density normalized to the nuclear matter density is shown 
in (3) for the lowest eigenmode. The time-like and transverse screening mass are seen to 
decrease midly with density. The longitudinal screening mass increases moderatly with 
baryon density. The mild dependence on the density for the SS model indicates that the 
vector mesons are weakly affected by the baryon density in this version of the SS model. As 
the inserted baryons are point like, at large N c their interaction is chiefly repulsive through 
w's as induced by D8-D8. The oo interactions with vectors and axials is mostly anomalous 
(through the WZ term) and therefore small as we ignored the WZ term. 



26 




6/J 

a 
"3 

<u 
o 
u 
o 
c/2 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



Longitudinal 
Transverse 



2 4 6 8 10 
nB / n 

(b) 



Figure 3: (a) Mass (b) Screening mass (Longitudinal mode: M' ] l (n j B)/m p (0), Transverse 
mode: M^(n B )/m p (0)) 



6.5.2 Pion decay constant 

The pion decay constant is identified from ((100), (105)) respectively, 

J n 



J dZK^A- 
AkMI 



The explicit dependence on the baryon density is shown in Fig. (4). Both the time-like and 




nB /n 



Figure 4: Pion decay constant 



space-like pion decay constant are found to increase with the baryon density. The increase 
is quadratic at small densities. Since the S-wave pion scattering length with baryons is 1/N C 
this explains the absence of a linear term. Moreover, for point-like external baryon sources 
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the pion- Axial- Vector coupling in matter at the origin of the pion decay constant involves 
two baryon sources and is repulsive. 

6.5.3 Vector Couplings and KSRF Relation 

The vector couplings are identified in (108) and (110). Their overall dependence on the 
baryon density is again mild as explained above. 



vlin couplings: 



g T 2 = 1 5 dZ 7tfe 

f dZ 

s _ 1 S dZ la 
9 V TP - cr > 7 i > 



(127) 



vvv couplings: 



T 
9 V 3 



9 V 3 



1 

K 



g s v ^^=JdZK-^{n x f^-\ 

dZR-^^fffA- 1 . (128) 



1 

K 




02468 10 02468 10 

n B / no n B / n 



(a) (b) 
Figure 5: (a) villi coupling (b) vvv coupling 
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6.5.4 KSRF relations 



In the matter rest frame Lorentz symmetry is no longer manifest. As a result, we expect a 
variety of KSFR relations depending on wether time-like or space-like parameters are used. 
Indeed, for instance the a-parameter at the origin of the KSFR relations can now take 4 
different forms depending on the time-like/space-like arrangement. Specifically 



Tl 

KSRF 



; 51 

KSRF 



4(^) 2 a 



T\2 



4(^) 2 (/; 



S\2 



(Mf) 



4(^) 2 (/J) 2 

4G40 2 (/, 5 ) 2 



,T2 
%SRF 
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'KSRF 



(129) 



2.2 
2 
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-S1 

a KSRF 






y 








-S2 
" 3 KSRF 
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4 6 
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Figure 6: Generalized a-parameter 



7 Conclusions 

We have considered a generalization of the chiral model proposed by Sakai and Sugimoto to 
finite baryon density. The baryon vertices in bulk are attached equally to the D8-D8 branes 
and correspond to S 4 in D8. They are treated as stable and point like in IR 3 and act as 
uniform sources of baryon density. Their point-like nature at large iV c and coupling A imply 
that their interactions as induced by D8-D8 is mostly repulsive through the exchanges of 
omega mesons. 

The bulk energy density grows quadratically with the baryon density before softening at 
asymptotic densities. The quadratic and repulsive growth is expected from the exchange of 
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omega mesons. The softening reflects on the fact that at asymptotic densities the repulsive 
baryons form an instable but regular array for fixed volume V. If V acting as a container 
is removed, the baryons fly away in this version of the SS model. We note that the energy 
density scales as N c since N c / y/a is of order 1 as expected from standard large N c arguments. 
The DBI action resums (partially) the strong NN-interactions while keeping the leading N c 
result unchanged. Since the instanton size is of order 1 / \/A we also note that the resummed 
contributions are of order A since the bulk instanton density \/\n B is of order A 2 (The 
additional \/\ here stems from the rescaling of z — ■> z/y/X in the delta-function source at 
z = 0). 

Using linear response theory, we have probed this dense baryonic system using pions, 
vectors and axials. The point like nature of the baryons with a size of order 1 / y/X and the 
large iV c nature as noted above, causes rather mild changes in the masses and couplings 
as a function of baryon density. In contrast, the pion decay constants are found to change 
appreciably. The quadratic increases at small baryon densities is mediated by omega's. The 
scalar S-wave pion-baryon scattering length is noted to vanish at large N c , causing f n to 
increase instead of decreasing at finite density. This behaviour is unphysical. 

The current approach needs to be improved in a number of ways to accomodate the 
baryon physics expected in the real world. First, the point-like nature of the sources need 
to be relaxed. This is possible by constructing the pertinent instanton vertex. Also, the 
point-like limit suggests that the DBI results quoted here are only indicative since higher 
derivative corrections to the DBI effective action are expected to contribute (see also [3, 4, 5] 
for further comments on this point)). Second, the Fermi motion of the sources need to 
be included. This can be achieved through a select quantization of the collective variables 
associated to the baryon vertex insertion. Some of these issues will be addressed in later work. 

Note added. After the completion of this work, we became aware of the recent work by O. 
Bergman, G. Lifschytz, and M. Lippert [15] who also address the SS model at finite baryon 
density. They have shown that a cusp configuration develops at finite density for generically 
separated D8-D8. This observation does not apply to the original SS embedding we discuss 
here. We also noticed the appearance of two relevant papers: [16] discusses the finite density 
problem in the holographic NJL model, and [17] discusses the effects of a finite size baryon 
charge distribution. 
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A The Effective Action Using the Vacuum Modes 

Let us only consider the lightest vacuum meson modes corresponding to the pion and p meson 
fields. This vacuum mode decomposition was studied in [3] at zero baryon density. Here we 
just add A as obtained in section 3 to the gauge field A M . Since the mode decomposition 
is complete, this approach should be complementary to the one discussed in the text. It is 
the same as the one we used in [8]. As we will show, the results are overall similar to the 
ones discussed in the main text regarding the density dependence. 
In the gauge A z = and £ = e ^ (54), A^ reads 

A^,z) = -iAo(z) +iv(OVi(2) . 

+ 0^n+ [^n 3 ]) Mz) + ^[rhAn] + o(n 4 ) , (130) 

where A is the background field, = We have set = for n > 2. The 

corresponding field strengths are 

2i ~ 1 ~ 

+ -([ ( 9 M n,^] + K, ( 9,n])^o + 7i [^n,^n](i-4^ 2 ) + o((n,^) 3 ) (m) 

J IT In 

F z , = -iAo+^jd^U+i^U^^ + v^ + OiU 4 ) (132) 

where A = ^ = % and 

0o = d g $ = ~ 0o in (37) (133) 



7 In this section the gauge field is treated as anti-Hermitian. Aq and II is Hermitian so i was introduced, 
while w M is anti-Hermitian. Note that we are working in a different gauge from Section 6. 
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Notice that A does not contribute to F^ and affect only F ZfJ/ . 
In order to compute the DBI action (39), 



R 3 9 U 

tr \l 1 — ^na'y—F^ - (2x0/)*-— F^F** + [F 3 ] + [F 4 ] + [F 5 ] , 

we need to know F^ U F^ U , F^ Z F^ Z , [F 3 ], [F 4 ], and [F 5 ], which have many complicated contri- 
butions. Again, we use the observations noted in the text to simplify. Thus 

[F A ] = ( W ) 4 ^ (fy* FtoF^FyF* + 0(K, y?) 4 ) • (134) 

Table (3) lists all relevant terms. We have introduced f^ v defined as 

Uu = dpv v - dyVp , (135) 

with fiu = 0,1,2,3 and i,j = 1,2,3. Table (3) should be understood in the integral and 
trace operation. We have omitted some terms vanishing in the operation and rearranged 
some terms by using the cyclicity of the trace. 

In terms of the entries in the RHS of the table, the action reads 



qDBI_ 
°D8-D8 



-T J d 4 xdzU 2 tr v 7 Po + Pi , (136) 



with 



P = 1-(2W) 2 ^ p = l-bK^d z A ) 2 , (137) 
4 Ukk 

Pi = (2™') 2 ^(a 2 + a 3 ) + (2 TO ') 2 ^(A+/3 2 ) 
2U 3 4 U KK 

+ (2»') 4 ^(7 2 + 73). (138) 
Again, P does not include meson fields and has carries the baryon density. Expanding the 
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= ^2 






= ^3 




- (A)) 2 




= A) 






--^((9 o n)0 o ^o + 2^ ^iA 


= /3i 

— Kl 








f/^i + fR n .^}Wi 


— t->2 
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l- 1 J 






— (2 






2/^[^,^](yio) 2 ^ + ^n](X) 2 ^(i - # 2 ) 












= 73 



Table 3: The relevant terms in evaluating DBI action up to third order in the fields (II, v). 
All entries are understood in the integral and trace operation. 



action by fluctuating the fields we have 



qDBI_ 
°D8-D8 



-fj 



(Txdz U 2 tr 



1 Pi 



1 Pi 
8 A /R? 



S 1 + S 2 + 0((U,v,) 4 ) , 



1 P[ 

i6 7n 



+ o((n,^) 4 ) 



(139) 



with 



S 1 = -T J d 4 xdz Uhr A' 1 , 



S 2 = -T / d A xdz U 2 ti 



-APi - -A 3 P 2 + — A 5 P 3 
2 1 8 1 16 1 



where we defined a modification factor A(Q) as 

1 1 



yf l-bK*(d z Ao) 2 



rat 



1 + ^- 5/3 - 



(140) 
(141) 



Notice that — Si is the grand potential discussed in section 3, and S 2 will be reduced to the 
action of mesons. To accomplish it we plug (138) into (141) and evaluate all z integrals and 
identify them as coefficients of each term in the remaining 4-D action. 
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dnpffiripnt 1 ^ 


T^pfinition 




T 


i fdZ^A 3 
i fdZK^A 


1 
1 






1 
1 


mi 
mi 


Ml K fdZK(d z ^ 1 ) 2 A :i 
M^ K J dZK(d z *i) 2 A 


m 2 p 
m 2 p 


T 

a 3 


jzJdZK-iUflA 
-± jdZK-HlA- 1 


4= • 0.446 
4 • 0.446 


T 
a ull2 

%n 2 


*4^/^^- 1/3 *i(l-4«)A 


v^4JV4 K - 1 - 584 



Table 4: The definitions of the coefficients in the action (142). At finite density, there are 
enhancing factors A, A 3 and a suppressing factor A -1 . 

Let us first check which terms we have and how they are affected by finite baryon density 
schematically. It can be read off from Table (3). At zero density we set A = and A = 1. 
Then a 2 ,a 3 ,p 2 survive. 0:2,^2 correspond to the free action of II and p, and a 3 is the 
couplings of vvv, villi interaction. At finite density all terms are enhanced by A, A 2 , or A 3 . 
Furthermore there are nontrivial modification. The free action part will be affected by 72 
and (3 2 . {f3\ itself does not contribute because the first term has odd parity in z and the 
second term is traceless.) The couplings are modified by 73. There are new interaction terms 
such as Vo(d fl v v — d u v^) 2 , VoV^, vod fl ILd fJ, IL, d ll{d^ll, v^ 1 }, which all vanish at zero density. 
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Coefficients 


Definition 


Q = 0(A = 1) 


T 




-4- -0.63 
-±- -0.63 



Table 5: The definitions of the coefficients in the action (142). At finite density, there are 
enhancing factors A, A 3 and a suppressing factor A -1 . 

Considering all these modification we get the final form of the meson action 



So = / d A x 



a£ 2 tr (9 n9°n) - 4 2 tr ($110*11) 
1 



+ a 1 v2 tTf 0i r + -a^trf ij f 
+ m 2T tr v v° + m^ 5 tr v iV % 



*3 



+ a 
+ 



> (2/ 0i [t; ,T; i ])+ aj.tr (/«[«*, t^"]) 
a^ n2 tr (2/ ,[«9°n, <9TI]) + a* n2 tr ^[0% ^n]) 



+ 



(142) 



where the coefficients of every term are defined in Table (4). At zero density all coefficients 
agree with those in [3]. There are three types of modification due to the baryon density: 
A, A -1 , and A 3 . A simply comes from \APi in (141). Since A is common to all coefficients, 
A 3 and A -1 can be understood as A ■ A 2 and A ■ A~ 2 . An enhancing factor A 2 is due to 
additional contribution from (3f and a suppressing factor A -1 is from 72,73, which explains 
the calculational similarities in all the results. 

Finally we want to mention without details, the character of the action in the gauge 
Am( z —> 00) — > instead of A z = 0. The action is the same as in (142) except for the 
interaction term v 1111 



-2^ n2 v [U, d°U] - 2g s vn2 Vi [u, dm] , 
where g^ n2 and g^ u2 are defined in Table(5). 



(143) 
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A.l Numerical results 



In this section we compute the coefficients in Table (4,5) numerically. Their physical mean- 
ings can be read off from the action (142). We use the same numerical inputs as detailed in 
section 3. 

A. 1.1 Pion Decay Constant 

The pion decay constant can be defined by the procedure of Section 6.3 with the vacuum 
mode function. 




(144) 



A. 1.2 Velocity 



The pion velocity: 




(145) 



The lowest mode velocity: 




(146) 



A. 1.3 Mass 




(147) 



A. 1.4 Screening mass 




(148) 
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